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ABSTRACT 

In the outer regions of the habitable zone, the risk of transitioning into a globally frozen 
"snowball" state poses a threat to the habitability of planets with the capacity to host 
water-based life. Here, we use a one-dimensional energy balance climate model (EBM) 
to examine how obliquity, spin rate, orbital eccentricity, and the fraction of the surface 
covered by ocean might influence the onset of such a snowball state. For an exoplanet, 
these parameters may be strikingly different from the values observed for Earth. Since, for 
constant semimajor axis, the annual mean stellar irradiation scales with (1 — e 2 ) -1 / 2 , one 
might expect the greatest habitable semimajor axis (for fixed atmospheric composition) to 
scale as (1 — e 2 ) -1 / 4 . We find that this standard simple ansatz provides a reasonable lower 
bound on the outer boundary of the habitable zone, but the influence of both obliquity and 
ocean fraction can be profound in the context of planets on eccentric orbits. For planets with 
eccentricity 0.5, for instance, our EBM suggests that the greatest habitable semimajor axis 
can vary by more than 0.8 AU (78%!) depending on obliquity, with higher obliquity worlds 
generally more stable against snowball transitions. One might also expect that the long 
winter at an eccentric planet's apoastron would render it more susceptible to global freezing. 
Our models suggest that this is not a significant risk for Earth-like planets around Sun-like 
stars, as considered here, since such planets are buffered by the thermal inertia provided by 
oceans covering at least 10% of their surface. Since planets on eccentric orbits spend much 
of their year particularly far from the star, such worlds might turn out to be especially good 
targets for direct observations with missions such as TPF-Darwin. Nevertheless, the extreme 
temperature variations achieved on highly eccentric exo-Earths raise questions about the 
adaptability of life to marginally or transiently habitable conditions. 

Subject headings: astrobiology - planetary systems -radiative transfer 



1. INTRODUCTION 

There are now more than 460 extrasolar planets 
known. 1 Selection effects favor the discovery of mas- 
sive giant planets orbiting very close to their parent 
stars, but advances in imaging and spectroscopic ca- 
pabilities, longer baselines for observation, and mis- 
sions such as CoRoT and Kepler should accelerate 
the rate of discovery of less massive planets in longer 
period orbits in the near future. In particular, mi- 
crolensing observations have already detected plan- 
ets less than ten times as massive as the Earth at 
distances as great as 2.6 AU (Beaulieu et al. 2006; 
Bennett et al. 2008), and have discovered a poten- 
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tial analog to our solar system (Gaudi et al. 2008) 
that might allow a habitable planet on a stable orbit 
(Malhotra & Minton 2008). Both CoRoT and Ke- 
pler promise to further increase the number of de- 
tected terrestrial exoplanets (Baglin 2003; Borucki 
et al. 2003, 2007; Borucki & for the Kepler Team 
2010). The observed secondary eclipse of HAT-P- 
7b by Kepler demonstrates that it should be capable 
of detecting transits of Earth-size planets (Borucki 
et al. 2009). Additionally, the CoRoT team recently 
announced the detection of the ~1.7R® exoplanet 
COROT-Exo-7b orbiting a KO star in the constel- 
lation Monoceros (Rouan et al. 2009; Leger et al. 
2009; Bouchy et al. 2009; Fressin et al. 2009) and the 
MEarth Project detected the transit of the 6.55 M e 
planet GJ 1214b (Charbonneau et al. 2009). Al- 
though these planets are in extremely close orbits 
(a = 0.017 AU, P = 0.85 days for COROT-Exo-7b; 
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a = 0.014 AU, P = 1.58 days for GJ 1214b) and are 
unlikely to be habitable, their discoveries represent 
a tremendous advance in planet detection capabil- 
ity. As the CoRoT, Kepler, and MEarth projects 
continue, even less massive terrestrial planets will 
surely be discovered in orbits with greater semima- 
jor axes. Once these potentially habitable terres- 
trial planets are discovered, researchers will be able 
to determine the radii, orbital semimajor axes, and 
masses of the planets, but current techniques are 
insufficient to constrain their obliquities and rota- 
tion rates (Valencia et al. 2006, 2007; Adams et al. 
2008). Although transit measurements might place 
constraints on the eccentricities of some transiting 
planets (Barnes 2007; Ford et al. 2008), radial ve- 
locity measurements of exact Earth analogs will be 
extremely challenging, and so the eccentricities of 
most such planets will remain undetermined in the 
near future. This paper attempts to quantify the 
effects of orbital parameters such as eccentricity 
on planetary habitability in order to prepare us to 
draw inferences about the habitability of yet-to-be- 
detected terrestrial planets even if their eccentrici- 
ties are not well constrained. 

However, current surveys of extrasolar planets 
indicate that the near-zero eccentricities seen in the 
Solar System are not necessarily typical and that 
many planets have significantly higher eccentrici- 
ties (Udry & Santos 2007). Of the exoplanets with 
measured eccentricities, ~40% are on more eccen- 
tric orbits than Pluto (e = 0.2488) and ~10% are on 
orbits with eccentricities >0.5. This suggests that 
current views of habitability that focus on direct 
Earth analogs in near-circular orbits might consider 
only a small subset of potentially habitable worlds. 
It therefore seems prudent to expand our study of 
habitability to encompass a wide range of orbital 
parameters. 

The study of planetary habitability began 
decades prior to the detection of exoplanets with 
the classic work of Dole (1964) and Hart (1979). In 
the last two decades, this topic has been revisited 
with increasing frequency, beginning with the work 
of Kasting et al. (1993), who found conservative lim- 
its for the Earth's liquid water habitable zone be- 
tween 0.95 AU and 1.37 AU. In recent years, various 
studies have applied the tools and techniques used 
to study the Earth's climate to simulations of plan- 
ets around other stars. Investigations by Williams 
& Pollard (2003), Williams & Kasting (1997) and 
Spiegel et al. (2009) have shown that, while varia- 
tions in the polar obliquity angle can alter the dis- 
tance from the star at which a planet becomes too 
cold to be habitable, planets with high obliquities 
are not necessarily less habitable than planets with 
low obliquities. Williams & Pollard (2002) also con- 
sidered the effect of eccentricity on habitability but 
only in the context of Earth twins at 1 AU. 

Spiegel et al. (2008, 2009) have analyzed the ef- 



fect of changes in semimajor axis, rotation rate, 
obliquity, and ocean coverage on the temperature 
of a generic terrestrial planet, but not variations 
in eccentricity. This paper attempts to fill in the 
void between the work of Spiegel et al. (2009) and 
Williams & Pollard (2002) by varying the eccentric- 
ity of model planets that are more diverse than the 
Earth- twin used by Williams & Pollard (2002). In 
addition, we also conduct a sensitivity study to de- 
termine which parameters have the strongest effect 
on planetary temperatures, so as to quantify the 
degree to which uncertainty in parameter measure- 
ment translates to uncertainty in climate. 

Although some planets may be habitable at all 
latitudes during their entire orbit, other planets, 
like the Earth, might be only partially habitable. 
These planets could therefore transition to a "snow- 
ball state" if small changes in insolation or atmo- 
spheric composition cause part or all of the planet 
to freeze. During a "snowball transition," the for- 
mation of snow or ice increases the albedo of the 
planet and the planet consequently becomes even 
colder. If the positive feedback loop between ice 
formation and increased albedo continues, the en- 
tire planet may become frozen and trapped in a 
snowball state. As discussed in Section 5.1, the ac- 
cumulation of greenhouse gases in the atmosphere 
while the planet is frozen may warm the surface suf- 
ficiently to allow the planet to eventually exit the 
snowball state. The transition from a snowball state 
to a partially habitable state is beyond the scope of 
this paper, but an investigation is pursued in Pier- 
rehumbert (2005) as well as in a companion paper 
(Spiegel et al. 2010). 

In this paper, based on a simple energy balance 
model (EBM) treatment, we determine that obliq- 
uity, eccentricity, and ocean fraction can together 
have a very strong influence on the orbital location 
of the snowball transition. For instance, for mod- 
els with an Earth-like atmosphere and eccentricity 
0.5 (which might not be extreme by extrasolar stan- 
dards), we find that the maximum habitable semi- 
major axis can extend to 1.90 AU or be as close 
to the star as 1.07 AU, depending on obliquity and 
ocean fraction. As a result, for e = 0.5, the stan- 
dard ansatz for the outer boundary of the habitable 
zone derived from considering only the annual mean 
flux (Barnes et al. 2008) can be off by more than 
78%. Altering the azimuthal obliquity (the degree 
of alignment of periastron with the solstices) does 
not have a significant effect on where the snowball 
transition occurs for low eccentricity planets, but 
the transition for high eccentric planets is pushed 
out significantly for azimuthal obliquity ~30° . More 
generally, planets in higher eccentricity orbits dis- 
play more latitudinal variation and seasonal vari- 
ation in habitability than planets in near-circular 
orbits. Because all of the models in this study in- 
corporate an Earth-like atmosphere, simultaneous 
variations of atmospheric composition or other fac- 
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tors in combination with the parameters considered 
in this study may produce a more complicated pic- 
ture of general planetary habitability. 

In Section 2 we discuss several important fac- 
tors that influence the Earth's climate on long time- 
scales. We explain the setup of our model in Section 
3 and discuss the validation of the model in Section 
4. In Section 5 we present our results. We then 
conclude in Section 6 and consider the implications 
of our findings on planetary habitability. 

2. GENERALIZED MILANKOVITCH CYCLES 

Milankovic (1941) realized that the long-term 
climate behavior of the Earth could in part be ex- 
plained by considering the combined effects of obliq- 
uity, orbital eccentricity, and precession. Each of 
these orbital elements changes on multiple, noncon- 
stant timescales known as Milankovitch cycles and 
the combination of these variations alters the cli- 
mate of the Earth by increasing or decreasing the 
solar insolation received by the Earth at a given lat- 
itude. For example, increasing the Earth's obliquity 
increases the annually averaged insolation received 
by the poles and decreases the annually averaged 
insolation received by the equator; both effects act 
to decrease the latitudinal temperature gradient." 
However, the Earth's obliquity is largely stabilized 
by the Moon and varies by only ±1.2° on timescales 
of -41 kyr (Laskar et al. 1993; Berger 1976, 1978). 
Much higher variations in obliquity are expected for 
planets without large moons: the obliquity of Mars 
varies between 14.9° and 35.5° (Ward 1974) and 
numerical models suggest that the Earth would ex- 
perience obliquity oscillations between 0° and 85° 
in the absence of the Moon (Laskar et al. 1993; 
Laskar & Robutel 1993). Even in the absence of 
large moons, however, the obliquity of a quickly 
rotating planet (< 8-hour day) would probably be 
self-stabilized by the fast rotation rate of the planet 
(Ward 1982; Laskar et al. 1993). 

In addition to obliquity, there is also a Mi- 
lankovitch cycle governing precession on shorter 
timescales of —19 kyr and —23 kyr (Berger 1976, 
1978). Over time, the slow shift of the direction 
of the Earth's rotation axis due to precession of 
both the spin axis and the orbital ellipse alters the 
position of solstices and equinoxes with respect to 
apoastron and periastron. In this paper, we con- 
sider this effect by varying the azimuthal obliquity 
angle a of our model planets (defined in Section 
3 as the angle between the position of the planet 
at periastron and the position of the planet at the 
northern winter solstice). Due to the ice albedo ef- 
fect, the hemisphere that is tilted toward the star 
at apoastron — which has a shorter winter (i.e., a 
shorter period of high albedo) — absorbs more inte- 
grated stellar energy per year than does the hemi- 
sphere that is tilted away at apoastron. Accord- 
ingly, the greatest temperature asymmetry between 
the northern and southern hemispheres is produced 



when periastron is aligned with a solstice (0 a = 0° 
or 6 a = 180°). Conversely, both hemispheres absorb 
equal annually averaged stellar irradiation when pe- 
riastron is aligned with an equinox (0 a = 90° or 
a = 270°). 

Finally, there is a Milankovitch cycle for ec- 
centricity. The Earth's orbital eccentricity is cur- 
rently nearly circular (e = 0.0167), but varies slowly 
up to —0.06 over long timescales of —100 and 
—400 kyr (Berger 1976, 1978). As shown in Figure 
1, the annual mean flux (F) scales as (1 — e 2 ) -1 / 2 
so that flux increases with increasing eccentricity 
for a given semimajor axis. Increasing the eccen- 
tricity, therefore, accentuates the ratio of the ir- 
radiating flux at periastron to that at apoastron 
(^max/^minOc((l + e)/(l-e)) 2 ), and slightly in- 
creases the annually averaged irradiation. The ratio 
^max/^min can be quite substantial for highly eccen- 
tric orbits, exceeding 10 2 for e > 0.82, which could, 
depending on a planet's thermal inertia and redistri- 
bution of energy, cause dramatic seasonal tempera- 
ture swings. We present climate models of planets 
with high orbital eccentricity in Section 5.1, and 
refer readers to a companion paper (Spiegel et al. 
2010) for a discussion of scenarios in which the ec- 
centricity of an Earth-like planet might be excited 
to large values. 

Over time, the combined changes in obliquity, 
precession, and eccentricity, together with nonlin- 
ear amplifications, can dramatically alter an exo- 
Earth's climate, leading alternately to periods of 
glaciation and of deglaciation as has been shown in 
the case of the Earth (Berger et al. 2005; Crucifix 
et al. 2006; Laskar et al. 1993; Loutre et al. 2004; 
Quinn et al. 1991). 

3. MODEL SETUP 

In this study, we investigate the temperature 
of a planet using the same one-dimensional time- 
dependent energy balance model introduced in 
Spiegel et al. (2008) and further explored in Spiegel 
et al. (2009). The model, which is similar to the 
more Earth-centric model used by Suarez & Held 
(1979), treats the meridional transport of heat as 
diffusion driven by the zonal mean temperature gra- 
dient: 

(i) 

This equation describes the evolution of the tem- 
perature T at location x = sin A, where A is the 
latitude, as a function of an effective heat capac- 
ity C, a diffusion coefficient D, and an albedo A. 
The net radiative energy flux in a latitude band is 
determined by the relationship between the energy 
received due to the diurnally averaged stellar flux 
S and the energy lost due to infrared emission /. 
One-dimensional EBMs such as this model provide 
a reasonable approximation of seasonal mean tern- 
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peratures for planets that rotate sufficiently quickly 
relative to their orbital frequency (Showman et al. 
2009). In all of our models we assume that the 
planet orbits a Sun-like star, so the stellar flux S 
is equivalent to that from a 1 M©, 1 L star. The 
effective heat capacities of the atmosphere over land 
(C/), over the wind-mixed surface layer of the ocean 
(C Q ) , and over ice (Ci) are the same as in Spiegel 
et al. (2008, 2009) and Williams & Kasting (1997) 
and are shown in Table 1. 

Previous work (Spiegel et al. 2008) explored 
three sets of infrared cooling radiation functions and 
albedo functions. Of the three models tested in that 
paper, Model 2 produced climates most similar to 
those on current Earth and is the one used in the 
current study: 



1 + 0.5925 (^k) 

where a is the Stefan-Boltzmann constant. The de- 
nominator of equation 2 would instead be 1 if the 
atmosphere had no opacity to outgoing infrared ra- 
diation; the functional form used for this denomina- 
tor represents the strength of the greenhouse effect, 
and is a reasonable approximation of the greenhouse 
for present-Earth conditions (Spiegel et al. 2008). 
In order to account for the higher reflectivity of ice 
and snow while using a simple functional form, we 
take the albedo to be constant and low (~0.28) at 
high temperatures, constant and high (~0.77) at 
low temperatures, and to vary smoothly in between: 

A[T] = 0.525 - 0.245 tanh f J . (3) 

The smooth hyperbolic tangent formulation is cho- 
sen to handle the phase transition from water to ice 
at 273 K in order to avoid the small ice-cap insta- 
bilities seen in models with a discontinuity in the 
albedo function at 273 K (Held et al. 1981). 

As explained in Spiegel et al. (2008), the 
fiducial diffusion coefficient follows the form 
of Williams & Kasting (1997) and is taken to be 
Am = 5.394 x 10 2 erg cm" 2 s" 1 K" 1 x (Q p /Q^)- 2 , 
where Q p and ^0 are the rotation rates of the 
exoplanet and the Earth, respectively. Its value 
increases with decreasing planetary angular spin 
frequency. We note that while the "thermal 
Rossby Number" scaling argument of Farrell (1990) 
supports this dependence of D on Q p , del Genio 
et al. (1993) and del Genio & Zhou (1996) find that 
the effective diffusivity of slowly rotating planets 
might not follow such a simple scaling relationship. 
The reader is directed to the recent review by 
Showman et al. (2009) for an in-depth discussion 
of the relationship between planetary rotation rate 
and effective diffusivity. 

The model is solved by relaxation on a grid of 
145 uniformly spaced latitude points using a time- 
implicit numerical scheme and an adaptive time- 
step, as described in Spiegel et al. (2008) and 



Hameury et al. (1998). We typically initialize the 
planet at the northern winter solstice, but for some 
runs we vary the azimuthal obliquity a to change 
the initial season. The initial temperature in all 
cases is uniform across the surface of the planet and 
is chosen to be at least 350 K to minimize the like- 
lihood that models evolve into ice-covered snowball 
Earths. Similar initial conditions were also used by 
Kasting et al. (1993) and Spiegel et al. (2008, 2009) 
to avoid "cold start" planets. As explained in the 
Appendix of Spiegel et al. (2008), as long as the ini- 
tial planet temperature is significantly warmer than 
273 K, our model relaxation studies indicate that 
within 130 years of model evolution (and sometimes 
far less), the final state of the planet is independent 
of the initial conditions. If the initial temperature 
is T < 273 K, however, the planet can quickly tran- 
sition to a snowball state due to water-ice albedo 
feedback. Our choice of initial conditions should 
therefore lead to more optimistic results for the lo- 
cation of the snowball transition. 

The purpose of this study is to examine the in- 
fluence of various orbital and planetary parameters 
on planetary habitability and to determine the sen- 
sitivity of a planet's habitability to changes in those 
parameters. Figure 2 portrays a schematic diagram 
indicating the relevant angles; here, we give a de- 
tailed description of parameters: 

1. Eccentricity e. We vary the eccentricity of our 
model planets from e = to e = 0.90. 

2. Polar Obliquity Angle 6 p . This is the angle be- 
tween the rotation axis of the planet and the 
normal to the plane of rotation. Because our 
model planets are symmetric, we restrict this 
angle to between 6 p = 0° and P — 90°. Val- 
ues between P = 90° and 6 P = 180° simply re- 
verse the designation of the identical northern 
and southern hemispheres. The wide range of 
polar obliquities is appropriate given the vari- 
ety of polar obliquities in our own solar system 
and the range of obliquities predicted by sim- 
ulations such as Agnor et al. (1999); Laskar 
et al. (1993); Laskar & Robutel (1993). In par- 
ticular, Kokubo & Ida (2007) suggest that po- 
lar obliquities near 90° may actually be more 
common than polar obliquities near 0°. 

3. Azimuthal Obliquity Angle 6 a . This is the an- 
gle between the projection of the rotation axis 
of the planet onto the plane of rotation and 
the line between the star and the periastron 
position of the planet. Variations in this angle 
change the initial season. The model is always 
initialized at periastron, and periastron coin- 
cides with northern winter for most models 
because most models have a = 0. Once the 
model reaches a periodic state, the average 
temperature of the planet will be greater at 
periastron than at apoastron. Consequently, 
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initializing the models at periastron is a rel- 
atively conservative choice because the aver- 
age planetary temperature will decrease as the 
planet approaches apoastron and the planet 
could enter a snowball state more quickly than 
a planet that was initialized with an aver- 
age global surface temperature of > 350 K at 
apoastron. 

4. Rotation Rate Q p . Because our parametriza- 
tion of the diffusion coefficient depends in- 
versely on the square of the rotation rate, in- 
creasing the rotation rate is equivalent to re- 
ducing the efficiency of latitudinal heat trans- 
port. In this study we consider planets with 
8-hour days (Q p = 3ft®, D = (1/9) Am), 24- 
hour days (Q p = ft®, D = Am) and 72-hour 
days (ft p = (l/3)ft®,Z) = 9AM). 

5. Semimajor Axis a. We begin each set of mod- 
els by using a simple scaling approximation 
as an ansatz about the location of the outer 
boundary of climatic habitability and then 
run a series of models with semimajor axes 
near that value to locate and refine the po- 
sition at which a model planet first becomes 
fully ice-covered year-round. 

6. Ocean Fraction f Q . This parameter deter- 
mines the ratio of land to ocean found on the 
surface of the planet. Since the atmosphere 
over the wind-mixed layer of the ocean has 
a much higher effective thermal inertia than 
the atmosphere over land, water world planets 
with high ocean fractions experience less dra- 
matic temperature variations than do desert 
worlds with low ocean fractions. Although the 
wind-mixed layer of the ocean varies in depth 
from a few meters to a few hundred meters 
or more (Hartmann 1994), we assume a depth 
of 50 m for the study. Increasing the depth 
of the wind-mixed layer would enhance the 
effective surface heat capacity and lengthen 
the timescale on which temperature changes 
occur. Theoretical simulations by Marotzke 
& Botzet (2007) indicate that the depth of 
the wind-mixed layer increases dramatically 
as the Earth freezes over, so our use of a 
constant 50-m wind-mixed layer may increase 
the likelihood that a planet will transition to 
a snowball state. Indeed, our most Earth- 
like model (a = 1 AU, e = 0, 6 P = 23.5°, 
6a = 0°, fo = 70%) transitions to a snow- 
ball state when the stellar luminosity has been 
reduced to 0.995L©; at 0.996L®, more than 
29% of the surface is covered by ice. This is 
comparable to the maximum stable ice cover 
of 30% reported by North (1975), but sig- 
nificantly below the value of 55% found by 
Voigt & Marotzke (2009). In this study, we 
present results from simulations with ocean 



coverage ranging between 10% and 90%, but 
most of our models incorporate an Earth-like 
70% ocean fraction (f Q = 0.7). Considering 
a range of ocean fractions is important be- 
cause simulations of planetary formation in- 
dicate a wide range of possible water contents 
(Raymond et al. 2004). In all cases, the land 
and ocean are uniformly distributed over the 
planet so that each latitude band has the same 
percentage of land and ocean coverage. Due 
to the lower specific heat capacity of land com- 
pared to water, altering the land distribution 
to produce polar continents could provide ad- 
ditional stability against snowball states. See 
Spiegel et al. (2009) for a detailed discussion 
of climate models of planets with nonuniform 
ocean coverage. 

4. MODEL VALIDATION 

Spiegel et al. (2008) confirmed that our EBM 
works reasonably well for the Earth and reproduces 
the Earth's climate to a degree of accuracy suf- 
ficient for investigations of exoplanet habitability. 
With the exception of the north/south asymmetry 
in the Earth's temperature profile at latitudes south 
of 60° S caused by Antarctica, the model agrees with 
the Earth's observed temperatures in 2004, as com- 
piled by the National Center for Environmental Pro- 
tection/National Center for Atmospheric Research 
(Kistler et al. 1999; Kalnay et al. 1996). 2 Spiegel 
et al. (2009) verified that the model predicts the sea- 
sonal variations in the Earth's radiative fluxes, by 
comparing the model results to data from NASA's 
Earth Radiation Budget Experiment (Barkstrom 
et al. 1990). The model results also agreed with 
those of Williams & Pollard (2003). 

As the current study considers eccentricity vari- 
ations as well as obliquity variations, we reexamine 
the model relaxation time to ensure that the model 
run time of 130 years used in Spiegel et al. (2008, 
2009) is still sufficient for the model to reach a re- 
laxed state under the forcing conditions explored 
here. As shown in the Appendix, model runs reach 
a stable oscillatory state within 100 years of model 
evolution. Thus, running the model for 130 years 
ensures that the resulting temperature profile rep- 
resents a relaxed state of the planet, independent of 
initial conditions. 

As a further check, we compare our results 
to those of Williams & Pollard (2002). In that 
work, a latitudinally resolved EBM and the three- 
dimensional climate code GENESIS 2 are used to 
model the climate of an Earth-like planet in an ec- 
centric orbit with a semimajor axis of 1 AU. The ge- 
ography, atmosphere, and obliquity angles of their 
model planet were identical to those of the Earth. 
Figure 3 is our version of their paper's Figure 2: 

2 Annually, zonally-averaged temperatures vary little from 
year-to-year. 
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we conducted model runs of planets with obliq- 
uity angles P = 23.5° and a = 0°, semimajor axis 
a = 1 AU, and eccentricities e = 0.1, 0.3, and 0.4. 
While our results do not strictly reproduce those 
of Williams & Pollard (2002), the general shapes of 
the temperature curves are similar and the tempera- 
tures generally agree to within 5 K. This agreement 
between our EBM and both the EBM and the gen- 
eral circulation model used by Williams & Pollard 
(2002) gives us further confidence in the suitability 
of our EBM for the habitability studies presented 
below. 

5. STUDY OF HABITABILITY 

We present a suite of models designed to probe 
the maximum semimajor axis at which a planet re- 
mains habitable before transitioning to a snowball 
state. We view habitability as a continuous, rather 
than a discrete, property and consider both tem- 
poral and regional habitability (Spiegel et al. 2008, 
2009). We also follow convention by adopting the 
freezing and boiling points of water under 1 bar of 
atmospheric pressure as the lower and upper bounds 
on habitable temperature. While boiling tempera- 
tures may seem extreme, there are several hyper- 
thermophiles on Earth that can grow at temper- 
atures above 373 K (e.g. Kashefi & Lovley 2003; 
Takai et al. 2008) so even our definition of habit- 
ability may be conservative and Earth-centric. Re- 
gardless, for the purpose of this paper, regions of 
a planet that are at temperatures between 273 K 
and 373 K are considered habitable while regions 
outside that temperature range are considered not 
habitable. We note that some regions of a planet's 
surface may be habitable even when the rest of the 
surface is not or that a planet may be habitable for 
only part of a year. Accordingly, we refer to both 
the temporal habitability fraction (the fraction of a 
year for which a given latitude band is habitable) 
and the regional habitability fraction (the fraction 
of the surface that is habitable at a given time). 
A detailed description of these terms is provided in 
Spiegel et al. (2008). 

When a planet becomes globally frozen year- 
round, we say that it has fallen into a "snowball" 
state. Here, we examine the maximum allowed 
semimajor axis that our models can withstand be- 
fore falling into a snowball state. Recall that 
the models in this paper do not include longterm 
geochemical feedback processes that would tend 
to stabilize a geophysically active planet's climate 
against such a snowball transition. As proposed 
by Walker et al. (1981), the decreased efficiency of 
the carbonate-silicate weathering cycle at low tem- 
peratures should cause greenhouse gases from vol- 
canic eruptions to accumulate in the atmosphere 
and gradually warm a planet out of a near-snowball 
state. Models by Kasting et al. (1993) that incor- 
porated this negative feedback loop showed that 
including the effects of the carbonate-silicate cy- 



cle extends the outer edge of the Sun's habitable 
zone to at least 1.37 AU. Since our model does 
not incorporate this feedback loop, our simulated 
planets may be more prone to snowball transi- 
tions than actual planets. Nevertheless, short-term 
changes in forcing may induce a snowball transi- 
tion in far less time than the ^million year pe- 
riod that would be required to accumulate enough 
CO2 in the atmosphere to restore temperate condi- 
tions. The value of our fixed- atmosphere models is 
to probe circumstances in which short-term (desta- 
bilizing) feedbacks might induce a snowball "phase- 
transition" that overwhelms longer-term (stabiliz- 
ing) feedbacks. 

5.1. Probing the Outer Limits of Habitability 

We probe the outer limits of habitability with a 
variety of diagnostic tests. We examine the effect of 
increasing semimajor axis while holding e constant 
(Figure 4) and the effect of increasing e at a con- 
stant semimajor axis (Figures 5, 6, and 7). We also 
explore the relative influence of rotation rate and ec- 
centricity for model planets at range of semimajor 
axes (Figure 8). Figures 6, 7, and 8 show dramati- 
cally increased outer boundaries of habitability for 
some eccentric models. In Figure 6, a planet with 
e = 0.9 and P = 0° has nonzero habitability out 
to 2.85 AU; in Figure 7, a planet with e = 0.7 and 
P = 23.5° is partially habitable to 1.215 AU; and 
in Figure 8, a planet with eccentricity of only 0.5 
and P = 90° (perhaps the most likely obliquity) is 
partially habitable out to 1.90 AU. These models all 
have Earth-composition atmosphere. 

Consider for instance the models in Figure 4 
for a planet with eccentricity e = 0.6, rotation rate 
Vtp = f^0, obliquity angles P = 23.5° and a = 0°, 
and ocean fraction = 70%. For a semimajor axis 
of 1.025 AU, the planet is completely habitable 
throughout the course of a year, but shows a strong 
temperature asymmetry due to the alignment of 
northern "winter" with periastron. Interestingly, 
the northern hemisphere actually decreases in tem- 
perature during northern "summer." This occurs 
because, even though the northern hemisphere is 
pointed toward the star, the planet is at apoastron 
and receives much less insolation than when at pe- 
riastron (see Figure 1). The southern hemisphere 
faces away from the star at apoastron and becomes 
much colder than the northern hemisphere because 
it receives even less insolation during its long winter. 

If the semimajor axis of the planet is increased to 
1.050 AU, the southern hemisphere becomes so cold 
during the winter that a permanent ice cover devel- 
ops over the southern pole in late southern spring. 
If the semimajor axis is increased to 1.075 AU, the 
ice forms earlier in the year because the planet cools 
faster and spends less time near the star at perias- 
tron. In addition, the ice that develops on a planet 
at a = 1.075 AU also extends farther northward be- 
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fore melting at periastron than the ice formed at 
smaller semimajor axes. Once the semimajor axis 
reaches 1.100 AU, however, the southern region of 
the planet becomes so cold during southern win- 
ter and spring that it cannot warm sufficiently at 
periastron and remains below freezing year-round. 
In addition, a seasonal northern polar cap develops 
during northern winter despite the proximity of the 
star. When the semimajor axis is further increased 
to 1.125 AU, the magnitude of the ice- water albedo 
effect is so strong that the entire planet transitions 
to a snowball state because more ice is formed dur- 
ing the southern winter at apoastron than can be 
melted during southern summer at periastron. 

As expected from the discussion in Section 2, 
the maximum habitable semimajor axis moves out- 
wards with increasing eccentricity. Figure 5 dis- 
plays the planetary temperature at each latitude 
over the course of one orbit for planets in orbits 
with a = 1 AU and increasing eccentricities. At 
low eccentricity neither pole receives enough inso- 
lation to raise the temperature above freezing dur- 
ing any part of the year, but, as the orbit becomes 
more eccentric, the net annual irradiation received 
by the planet increases (see Figure 1) and the hab- 
itable region of the planet expands. The north- 
ern pole becomes habitable at a lower eccentricity 
than the southern pole because the southern hemi- 
sphere faces away from the star at apoastron and 
therefore absorbs less annually averaged flux than 
the northern hemisphere. The southern hemisphere 
also experiences a longer, colder winter than the 
northern hemisphere, which faces toward the star 
at apoastron and away from the star at periastron. 
The model planet's polar climate differs from that 
of the Earth because of a variety of simplifications 
that the model has compared to the Earth, includ- 
ing uniform distribution of continents and ocean, 
the azimuthal obliquity of 0° (in comparison to the 
Earth's value of 13°), the various climate feedbacks 
that our model does not include, and our treatment 
of heat redistribution as a purely diffusive process. 

A visual example of the outward movement of 
the maximum habitable semimajor axis with in- 
creasing eccentricity is provided in the left panel of 
Figure 6, which displays the temporal habit ability 
fraction of a series of planets with polar obliquity 
P = 0° and a range of eccentricities and semima- 
jor axes. As the eccentricity of the planet's orbit is 
increased, the maximum habitable semimajor axis 
also increases for each latitude band of the planet. 
The most extreme example is shown in the last 
row of Figure 6 for planets with e = 0.9. At such 
high eccentricities, the range of semimajor axes for 
habitable planets lies past the maximum habitable 
semimajor axis for planets in low-eccentricity orbits 
(e < 0.10). In addition, the variation of habitability 
with latitude is much more noticeable for the case of 
e = 0.9 planets than for the lower eccentricity plan- 
ets. At the equator, this model maintains partial 



habitability out to 2.85 AU. While this is an inter- 
esting suggestion that such highly eccentric planets 
could have dramatically expanded habitable zones, 
this result should be viewed cautiously, since such 
extremely eccentric planets might not be reasonably 
modeled by a simple EBM. 

Seasonal variations in habitability are much 
more pronounced on planets in highly eccentric or- 
bits. As shown in the right panel of Figure 6, the 
habitability of a planet in a near-circular orbit is 
nearly constant year-round, but the habitability of 
planets with e > 0.30 depends on the season. Be- 
cause the models were initialized at periastron, the 
planet receives much more stellar irradiation dur- 
ing the first half of the year (time < 0.5 years) than 
during second half of the year (time > 0.5 years). 
During the long winter near apoastron, ice accumu- 
lates at the poles and decreases the habitability frac- 
tion of the planet, especially at increased values of 
the semimajor axis. Near periastron, however, the 
seasonal ice cover melts and the regional habitabil- 
ity fraction is increased. Accordingly, the regional 
habitability fraction for planets on highly eccentric 
orbits depends strongly on both the semimajor axis 
and the time of the year. 

Since the model planets analyzed in Figure 6 
have P = 0° and a = 0°, their temporal habitabil- 
ity fraction is symmetric with respect to the equa- 
tor. For planets with non-zero p , the northern 
and southern hemispheres display different tempo- 
ral habit abilities. Figure 7 shows the temporal and 
regional habitability fractions for a set of model 
planets with P = 23.5° and a range of eccentric- 
ities. As displayed in the left panel, the northern 
hemisphere of such planets is more habitable than 
the southern hemisphere at larger semimajor axes. 
This effect is most pronounced for planets in highly 
eccentric orbits. 

The right panel of Figure 7 demonstrates that 
planets with eccentricities between 0.4 and 0.6 ex- 
perience sharp decreases in regional habitability at 
the snowball transition. Just inside the maximum 
habitable semimajor axis, at least 60% of the sur- 
face of a planet with a 24-hr day or a 72-hr day is 
habitable for the majority of the year. Less than 
0.005 AU beyond this distance, however, the en- 
tire planet becomes completely non-habitable. At 
higher eccentricities, the snowball transition is much 
more gradual. For eccentricity 0.7, a small region 
of the planet remains transiently habitable for semi- 
major axes between 1.2 AU and 1.225 AU regardless 
of rotation rate. Consequently, the regional habit- 
ability plot reveals a sharp cut-off in regional hab- 
itability at the snowball transition for planets with 
moderate eccentricities (0.4 < e < 0.6) but a long 
"tail" of decreasing partial habitability for slightly 
higher eccentricity (e = 0.7). 

Figure 8 shows the semimajor axis correspond- 
ing to the snowball transition for planets with 8-, 
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24-, and 72-hour days. The semimajor axis indi- 
cated is the largest semimajor axis for which at 
least part of the planet is habitable at some point 
in the year. For reference, the black curve plots 
the ratio of the annual mean flux received on an 
orbit with the eccentricity shown along the x-axis 
to the flux received on a circular orbit at 1 AU. 
As displayed in the figure, the assumption that the 
location of the habitable zone scales with the orbit- 
averaged flux is justified for low to moderate eccen- 
tricity orbits (e < 0.65), but for highly eccentric 
orbits, the habitable zone extends to much greater 
distances than that simple scaling relationship pre- 
dicts. For e = 0.75, for example, model planets 
with p = 23.5°, 9 a = 0, and f Q = 0.7 are hab- 
itable out to a ~1.39 AU even though the scaling 
relationship would predict an outermost habitable 
semimajor axis of only ^1.23 AU. Despite the sim- 
plicity of this scaling relation, our ansatz is typically 
within the identified snowball transition region for 
a planet with a 24-hour day. 

Intriguingly, in the context of our models, the re- 
lationship between rotation rate and the position of 
the snowball transition does not seem to be mono- 
tonic. In a circular orbit, a planet with an 8-hr 
day, P = 23.5°, a = 0, and f Q = 0.7 under- 
goes a snowball transition at a ~0.95 AU, but more 
slowly rotating planets (24-hr days or 72-hr days) 
are habitable out to a ~1 AU. Consequently, more 
quickly rotating planets in circular orbits freeze over 
at shorter distances than more slowly rotating plan- 
ets. However, if the eccentricity of the orbit is 
raised to e = 0.2, planets with Earth-like 24- hr days 
are habitable at greater semimajor axes (a ~1.013) 
than planets with either 8-hr days (a ~ 1.003) or 72- 
hr days (a ~0.993). These differences are smaller 
than the difference in the position of the snowball 
transition for quickly rotating (8-hr days) and less 
rapidly rotating planets (24-hr or 72-hr days) in cir- 
cular orbits, but the order of the distances of the 
snowball transitions is unexpected. As discussed by 
Spiegel et al. (2008), this suggests that there may 
be a trade-off in keeping the equator warm by rotat- 
ing sufficiently quickly that not all of the heat can 
diffuse to the poles and by rotating slowly enough 
to allow enough heat to diffuse to the high latitudes 
to prevent the formation of large-scale polar ice cov- 
erage that could cool the entire planet through the 
ice- water albedo effect. Exploring a denser range 
of rotation rates across a variety of eccentricities 
might help elucidate the relationship between rota- 
tion rate, eccentricity, and the position of the snow- 
ball transition in low eccentricity orbits. For more 
eccentric orbits, however, rotation rate (or at least 
the rotation rates studied here) does not appear to 
have a significant influence on the position of the 
snowball transition. As shown in Figure 8, for or- 
bits with e = 0.5, the onset of the snowball transi- 
tion occurs at the same semimajor axis for planets 
with 8-hr, 24-hr, and 72-hr days. Additionally, the 



position of the snowball transition is nearly identi- 
cal for planets with 24-hr and 72-hr days in orbits 
with e> 0.35. 

5.2. Sensitivity Study 

Previous models of habitability (Spiegel et al. 
2008, 2009; Williams & Pollard 2002, 2003; Williams 
& Kasting 1997) have investigated a variety of test 
planets, but there is still a large region of parame- 
ter space unexplored. The sheer number of factors 
influencing climatic habitability means that thou- 
sands of model runs would be required to fully 
explore the contours of the snowball transition in 
the multi-dimensional space of parameters describ- 
ing the star, the orbit, the planetary spin (rate 
and obliquity), and properties of the planet's atmo- 
sphere and surface. Instead, we present the results 
of a sensitivity study to determine which parame- 
ters have the strongest effects on habitability. We 
find that increasing the polar obliquity increases the 
semimajor axis of the snowball transition for both 
low and high eccentricity planets, but that the ef- 
fects of changes in azimuthal obliquity or ocean cov- 
erage depend on eccentricity. 

We consider two model planets, one with eccen- 
tricity 0.2 and the other with eccentricity 0.5. Both 
model planets have Earth-like obliquity (0 P = 23.5°, 
a = 0° compared to P(B = 23.5°, a(B = 13° for the 
Earth) and uniform continent distributions with an 
Earth-like 70% ocean fraction. We first determine 
the maximum habitable semimajor axis for both 
planets by conducting preliminary model runs on 
a fine grid in semimajor axis (0.005 AU spacing). 
Then, we systematically alter each parameter ei- 
ther individually or in combination to determine the 
maximum semimajor axis for a habitable planet as 
a function of slight deviations in each input param- 
eter. Our results for the planet with e = 0.2 are 
shown in Table 2 and our results for the planet with 
e = 0.5 are shown in Table 3. 

For both cases increasing the polar obliquity al- 
lows the planet to remain habitable at greater semi- 
major axes. When P < 57°, the equator receives 
more annually averaged insolation than the poles 
and the direction of heatflow is poleward, but at 
higher polar obliquities, the poles receive more inso- 
lation than the equator and the direction of heatflow 
is reversed. As a result, while planets with P = 0° 
develop permanent ice coverage at both poles and 
transition to a snowball state when the downward 
extension of the ice reaches the equator, planets 
with Op = 90° are warmest at their poles and heat 
is transported from the poles toward the equator. 
Planets with 90° polar obliquity therefore develop 
small, seasonal ice coverage near the equator and 
remain habitable for an additional 0.04 AU (for 
e = 0.2) or for an additional 0.838 AU (for e = 0.5). 
This later case is the planet that is habitable at 
1.90 AU featured in Figure 8, which is transiently 
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habitable at the south pole during the brief, intense 
periastron summer. 

Reducing the ocean fraction from 70% ocean to 
10% ocean decreases the maximum habitable semi- 
major axis for planets with high eccentricity, but 
has a more complicated effect on planets with low 
eccentricity depending on polar obliquity. Con- 
sider, for example, a desert planet with 6 p = 90°, 
a — 1.255 A, and e = 0.2. The desert planet expe- 
riences tremendous temperature oscillations during 
the course of the year, but the southern pole be- 
comes transiently habitable during northern winter 
at periastron. The southern pole then freezes dur- 
ing the long winter and reaches cold temperatures 
(T ~150 K) before thawing again at apoastron. As 
shown in Figure 9, the southern pole is the only 
part of the planet that is ever within the liquid wa- 
ter limits of habitability and experiences the most 
extreme temperature variations on the planet. How- 
ever, properly modeling this planet would require 
taking into account the latent heat of freezing and 
melting. 

Variations in azimuthal obliquity shift the po- 
sitions of the solstices with respect to periastron 
and apoastron. As a result, the azimuthal obliquity 
can influence the orbital distance at which a planet 
transitions to a snowball state. At low eccentric- 
ity, there is not much effect. At higher eccentricity, 
6 a ~ 30° leads to a larger habitable semimajor axis 
than 9 a = 0° or a = 90°. Why might this occur? 
At low azimuthal obliquity, the contrast between 
the annually averaged flux received by the south- 
ern and northern hemispheres is conducive to ice 
growth. Conversely, for a ~ 90° the planet experi- 
ences milder winters, but the less intense summers 
reduce seasonal melting. A sweet spot occurs near 
6 a = 30°, where the enhanced irradiation experi- 
enced during summer roughly compensates for the 
colder winters. As shown in Figure 1, the ratio of 
the flux received at periastron to the flux received at 
apoastron increases with increasing eccentricity so 
that variations in azimuthal obliquity have a more 
pronounced effect on planets in more eccentric or- 
bits. 

6. SUMMARY AND CONCLUSION 

We presented the results of a series of idealized 
energy balance model runs to determine the habit- 
ability of planets with a range of eccentricities. As 
shown in Section 5.1, planets in orbits with a given 
semimajor axis can remain habitable for a range of 
eccentricities. This suggests that if a planet were 
to experience eccentricity perturbations caused by 
a giant planet companion, as considered in the com- 
panion paper (Spiegel et al. 2010), it would not nec- 
essarily transition to a snowball state. Instead, our 
study suggests that many perturbed planets could 
remain fully or partially habitable even at slightly 
different eccentricities. Intriguingly, a previously 
frozen planet might thaw if perturbed to a higher 



eccentricity. Incorporating the latent heat involved 
in the melting of a frozen world is a complex prob- 
lem that we do not address in this paper, but we 
propose a solution in Spiegel et al. (2010). 

Throughout this study, we have used the con- 
ventional liquid water definition of habitability, but 
many extremophiles are capable of surviving outside 
of that temperature range (Carpenter et al. 2000; 
Kashefi & Lovley 2003; Takai et al. 2008; Junge 
et al. 2004). While it remains to be seen whether life 
could originate at temperatures well below 273 K or 
well above 373 K, we should remain cautious about 
making any assumptions about the limitations of 
microbes. Recent advances in biology continue to 
demonstrate that lifeforms on Earth are far more 
inventive and ubiquitous than we ever would have 
expected. 

In this study we have considered the effects of 
eccentricity on the semimajor axis corresponding 
to the snowball transition for pseudo-Earth planets 
and found, as has been seen by Williams & Pollard 
(2002), that increasing eccentricity can increase the 
allowed semimajor axis for habitable planets. In ad- 
dition to increasing the maximum semimajor axis at 
which the planet can remain habitable, increases in 
eccentricity also enhance regional and seasonal vari- 
ability in planetary temperatures and lead to a more 
gradual transition from habitable to non-habitable 
planets with increasing semimajor axis. 

We have also conducted a sensitivity study to de- 
termine which orbital parameters are the most influ- 
ential on the location of the snowball transition. Al- 
though changing the obliquity of a low-eccentricity 
planet can alter the maximum habitable semimajor 
axis by a hundredth of an AU, we find that reduc- 
ing the ocean fraction has the strongest effect on 
the maximum habitable semimajor axis of a low- 
eccentricity planet because of the important role 
of the ocean's thermal inertia in mediating climate 
variations. Combining decreases in ocean coverage 
with increases in polar obliquity can further extend 
the position of the snowball transition, but changes 
in azimuthal obliquity have a negligible effect on the 
habitability of low eccentricity planets. 

The situation for higher-eccentricity planets is 
more complicated, but variations in polar obliquity 
seem to have the most powerful effect on the posi- 
tion of the maximum habitable semimajor axis and 
can increase the semimajor axis corresponding to 
the snowball transition by ~0.8 AU. Changes in 
azimuthal obliquity are also significant for highly 
eccentric planets because of the uneven distribu- 
tion of flux throughout the orbit. Increasing the 
azimuthal obliquity to 90° actually decreases the 
semimajor axis of the snowball transition, but there 
is a sweet spot near ~30° where the increase in 
the azimuthal obliquity extends the position of the 
snowball transition by ~0.175 AU. Our simulations 
indicate that changes in the effective thermal dif- 
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fusivity by roughly an order of magnitude in either 
direction (motivated by the suggestion that thermal 
diffusivity might depend strongly on rotation rate) 
have little influence on habitability for planets with 
moderate eccentricity (0.35 < e < 0.65), but planets 
with high or low eccentricity display a complex de- 
pendence of the position of the snowball transition 
on diffusivity. 

Finally, this study suggests that planets in 
moderately- or highly-eccentric orbits (e > 0.5) may 
be habitable to much larger semimajor axes than 
would result from simply scaling the semimajor axis 
to match the flux received in a circular orbit. In 
particular, a model with e = 0.9 and 9 P = 0° is 
habitable to 2.85 AU and a model with e = 0.5 and 
P = 90° is habitable to 1.90 AU, both with Earth- 
composition atmospheres. Although these numbers 
are surprisingly large for a fixed-composition atmo- 
sphere, the fact that EBMs tend to be more sus- 
ceptible to global freezing than are actual plan- 
ets gives us some confidence that our models are 
not prone to overestimating the outer boundary of 
habitability. The partially habitable model with 
e = 0.5 and a = 1.90 AU has apoastron separa- 
tion of 2.85 AU, which raises the intriguing possibil- 
ity that some moderately-to-highly eccentric plan- 
ets in the outer reaches of a habitable zone might 
be significantly easier to observe with direct imag- 
ing platforms such as TPF-Darwin (Kaltenegger & 
Fridlund 2005; Heap et al. 2008) than similar plan- 
ets on circular orbits. Nevertheless, this question 



cannot be properly evaluated without a model that 
appropriately accounts for both the latent heat of 
melting/freezing and the atmospheric changes (in 
composition, cloud-cover, etc.) that occur with such 
strong changes in stellar irradiation over the an- 
nual cycle. Ideally, future studies would combine 
study of changes in polar obliquity, azimuthal obliq- 
uity, rotation rate, and eccentricity with variations 
in other parameters such as continent distribution, 
atmospheric composition, and azimuthal obliquity 
to further explore the multi-dimensional contours 
of the snowball transition. 
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APPENDIX 

A detailed model relaxation study was completed in Spiegel et al. (2008), but because the present work 
focuses on variations in eccentricity, we repeat the test to ensure that 130 years of model evolution remains 
sufficient for the planets in this study. Figure A shows the mean global temperature of several model planets 
smoothed with a 1-year boxcar filter. Two of the model planets have rotation rate Q p = ^0 and the third 
model has rotation rate Q p = (1/3)^0. In this case, the more slowly rotating planet reaches a steady climate 
state after a shorter model runtime than the planets with 24-hr days. However, as shown in detail in the 
right panel, all of the models achieve a steady climate state within 100 years and then show maximum 
temperature deviations of less than 0.04 K. Specifically, the model planet with Vt p = e = 0.15, and 
a = 1.005 AU reaches a steady climate state after ~17 years and the model planet with Vt p — fJ e , e = 0.75, 
and a — 1.125 AU reaches a steady climate state after ~90 years. The model planet with Vt p = (l/3)O e , 
e = 0.15, and a = 0.950 AU achieves a steady climate state after only ~5 years of model evolution. The 
results of this study indicate that integrating for 130 years is sufficient for this purpose and that the time 
required to reach a steady climate state does not increase monotonically with increasing eccentricity. 
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Fig. A. — Model relaxation tests. Left: Variations in the average global temperature over the 130-year model evolution 
for model planets with polar obliquity P = 23.5° and azimuthal obliquity a = 0° and a variety of eccentricities, rotation 
rates, and semimajor axes. The average global temperature is smoothed by a 1-year boxcar filter for the entire 130-year model 
evolution of a planet. The slowly rotating planet reaches a steady climate state after <10 years, but as shown in the right panel, 
the planets with 24-hr days require longer model runs to reach a steady climate state. Right: A more detailed plot of the 
average global temperature for the model planets with rotation rate Q p = f2® highlighting the transition to a steady climate 
state. The planet with eccentricity e = 0.15 cools to T ~290.3 K after ~90 years, but the planet with eccentricity e = 0.75 
reaches a steady climate state with T = 314.7 K after only ~17 years. As before, the global mean temperature is smoothed by 
a 1-year boxcar filter. 

TABLE 1 

Model Values for the Surface Heat Capacity C 



Surface Type 


Effective Heat Capacity (erg cm 2 K L ) 


Land 


C t = 5.25 x 10 y 


Ocean 


Co = 40Q 


Ice (263 K < T < 273 K) 


d = 9.2C Z 


Ice (T < 263 K) 


Ci = 2.0d 



TABLE 2 

Sensitivity Study for a Planet with Eccentricity e = 0.2 



Variation in Model Outer Edge Change from 

Parameters of HZ (AU) Fiducial Model (AU) 





Fiducial Model (Earth with e = 0.2) a 
= 23.5° a = 30° Ocean Fraction = 




1.010-1.015 


N/A 


6p 


0.7 


1.010-1.015 


0.000 


Op 


= 23.5° 


a = 90° 


Ocean Fraction = 


0.7 


1.010-1.015 


0.000 


Op 


= 0° 


a = 0° 


Ocean Fraction = 


0.7 


1.015-1.020 


0.005 


Op 


= 15° 


a = 0° 


Ocean Fraction = 


0.7 


1.010-1.015 


0.000 


Op 


= 30° 


a = 0° 


Ocean Fraction = 


0.7 


1.010-1.015 


0.000 


Op 


= 60° 


a = 0° 


Ocean Fraction = 


0.7 


1.030-1.035 


0.020 


Op 


= 90° 


a = 0° 


Ocean Fraction = 


0.7 


1.050-1.055 


0.040 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.9 


1.010-1.015 


0.000 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.5 


1.010-1.015 


0.000 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.1 


0.990-0.995 


-0.020 


Op 


= 0° 


a = 0° 


Ocean Fraction = 


0.1 


0.995-1.000 


-0.015 


Op 


= 90° 


a = 0° 


Ocean Fraction = 


0.1 


1.255-1.260 


0.245 



a The parameters for this planet are f Q = 0.70, P = 23.5, a = 0, and fl p = 
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TABLE 3 

Sensitivity Study for a Planet with Eccentricity e = 0.5 







Variation in 




Outer Edge 


Change from 






Model Parameters 




of HZ (AU) 


Fiducial Model (AU) 




Fiducial Model (Earth with e = 0.5) b 
= 23.5° a = 30° Ocean Fraction = 




1.070-1.075 


N/A 


6p 


0.7 


1.245-1.250 


0.175 


Op 


= 23.5° 


a = 90° 


Ocean Fraction = 


0.7 


1.065-1.070 


-0.005 


Op 


= 0° 


a = 0° 


Ocean Fraction = 


0.7 


1.070-1.075 


0.000 


Op 


= 15° 


a = 0° 


Ocean Fraction = 


0.7 


1.070-1.075 


0.000 


Op 


= 30° 


a = 0° 


Ocean Fraction = 


0.7 


1.365-1.370 


0.295 


Op 


= 60° 


a = 0° 


Ocean Fraction = 


0.7 


1.765-1.785 


0.703 


Op 


= 90° 


a = 0° 


Ocean Fraction = 


0.7 


1.900-1.920 


0.838 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.9 


1.075-1.080 


0.005 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.3 


1.090-1.095 


0.020 


Op 


= 23.5° 


a = 0° 


Ocean Fraction = 


0.1 


1.255-1.260 


0.185 


Op 


= 0° 


a = 0° 


Ocean Fraction = 


0.1 


1.070-1.075 


0.000 


Op 


= 90° 


a = 0° 


Ocean Fraction = 


0.1 


1.415-1.420 


0.345 



b The parameters for this planet are f Q = 0.70, P = 23.5, a = 0, and Q p = Qq. 
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Fig. 1. — Comparison of the mean flux (blue) and the ratio of periastron flux to apoastron flux (red) for planets in eccentric 
orbits to those of planets in circular orbits. The top plot is logarithmic and extends from e = to e = 1 while the bottom 
plot is linear to highlight the change in flux between e = and e = 0.5. For reference, the green line shows the distance of the 
planet at apoastron relative to the distance of the planet at periastron. Although the average flux changes by less than a factor 
of two even at very high eccentricity, the change in flux over the course of the orbit is several orders of magnitude larger for 
planets at high eccentricity. Consequently, regions of those planets may be both well below 273 K and well above 373 K during 
the course of a year. Provided the planets do not freeze over completely during the long winter and enter a snowball state, this 
raises the question of whether life would be able to withstand such extreme temperature variations. 
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Fig. 2. — Top view of a planet in an elliptical orbit with semimajor axis a and eccentricity e. As explained in Section 3, the 
azimuthal obliquity a is zero if periastron is aligned with the projection of the planet's spin angular momentum vector onto 
the orbital plane. In that case, northern winter coincides with periastron. 
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Fig. 3. — Our version of Figure 2 of Williams Sz Pollard (2002). Their model used a test planet with the geography of 
the Earth while our model used a test planet with a uniform continent distribution but the same 70% ocean fraction as the 
Earth. The slight differences between their results and ours might be due to the presence of an asymmetric southern continent 
(Antarctica) in their model but not ours, and other model differences between their three-dimensional GCM and our simpler 
one-dimensional EBM. 



18 




129.25 129.5 129.75 130 129.25 129.5 129.75 130 

Last 1 Year Last 1 Year 

Fig. 4. — Temperature as a function of latitude and time of year for a set of planets with e = 0.6, P = 23.5°, a = 0°, 
f Q = 70% and rotation rate Q p = Q§ with a range of semimajor axes. The color indicates temperature as shown in the colorbar. 
Areas of the planet shown in white are either below 273 K or above 373 K. The planet is completely habitable at 1.025 AU, but 
the temperature near the south pole drops below freezing during the southern winter at 1.050 AU. The southern ice coverage 
grows for semimajor axes > 1.075 AU and is accompanied by northern ice coverage during northern winter at 1.100 AU. Beyond 
1.125 AU, the planet is completely frozen. 
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129.25 129.5 129.75 130 129.25 129.5 129.75 130 

Last 1 Year Last 1 Year 

Fig. 5. — Temperature as a function of latitude and time of year for a set of planets with P = 23.5°, a = 0° , f Q = 70% and 
rotation rate Q p = All of the planets are in an orbit with semimajor axis a = 1 AU, but the eccentricity varies from e = 
(top left) to e = 0.15 (bottom right). At near-zero eccentricity, only the mid-latitudes are habitable, but at high eccentricities, 
the polar regions receive greater mean isolation over the course of the year and the regional habitability fraction increases. 
In addition, note that the temperature asymmetry between the northern and southern hemispheres due to the alignment of 
periastron with northern winter is more pronounced at higher eccentricities because northern winter occurs much closer to the 
star than southern winter. 



20 



e = 0.01 



e = 0.10 



e = 0.30 



e = 0.50 



e = 0.90 



0.6 0.85 1.1 1.35 1.6 1.85 2.1 2.35 2.6 2.85 
a (All) 



yu 






i 












CO 


45 


« 






CD 
>" 





1 




0.9 


o 








CD 


-45 


5 






Tim 


-90 




- 


0.8 












CO 


45 


CD 
"O 






CD 

> 







-- 


0.7 


o 


to 




CD 


-45 


i 






Tim 


-90 






0.6 












CO 


45 


CD 

-Q 






CD 

>- 







-- 


0.5 


o 




"co 




CD 


-45 


i 






Tim 


-90 






0.4 




45 








CO 


CD 






CD 

>- 







-- 


0.3 


o 








CD 


-45 








Tim 


-90 


: 




0.2 




45 




i 




CO 








CD 

>- 









0.1 


O 






CD 


-45 








Tim 


-90 












e = 0.01 



0.251 



0.51 
0.251 



e = 0.10 



e = 0.30 



e = 0.50 



^^^^^^e = 0.90 



0.6 0.85 1.1 



1.35 1.6 1.85 2.1 
a (All) 



2.35 2.6 2.85 



Fig. 6. — Left: Temporal habitability fraction (the fraction of a year for which a given latitude band is habitable) for planets 
with 6 P = 0°, 9 a = 0°, ocean fraction = 70%, and Q p = fi®. The eccentricity of each set of planets is indicated and increases 
downward from e = 0.01 to e = 0.90. Note that the maximum habitable semimajor axis increases for greater eccentricities. 
Regions shown in red indicate that the region is habitable for the entire year while regions shown in blue indicate that the 
region is non-habitable for the entire year. Colors between red and blue signify that the region is habitable for part of the year. 
Right: Regional habitability fraction (the fraction of the surface that is habitable at a given time) for the same planets as in 
the left panel. Note the increased seasonal dependence at higher eccentricities due to the extreme variations in flux over the 
orbit. Times shown in red indicate that the entire planet is habitable at that time while times shown in blue indicate that the 
entire planet is non-habitable at that time. Colors between red and blue signify that a fraction of the planet is habitable at 
that time. 



21 
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a (AU) a (AU) 

Fig. 7. — Same as Figure 6 but for planets with P = 23.5°. Eccentricity increases downward from e = 0.4 for the top 
plots to e = 0.7 for the bottom plots. Left: Temporal habitability fraction as a function of semimajor axis. Right: Regional 
habitability fraction as a function of semimajor axis. Note that the light shades of blue at a = 1.2 AU in the plots corresponding 
to planets with eccentricity e = 0.7 means that a small fraction of the planet is habitable for a brief period in during southern 
fall. Comparison with the left panel reveals that the area of the planet that is is transiently habitable during southern fall is 
near — 30° S, which is near the substellar point at periastron. 
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Fig. 8. — Effect of rotation rate on the position of the snowball transition. Each line displays the approximate location 
of the snowball transition for planets with 72-hr days (red), 24-hr days (green), and 8-hr days (blue). For comparison, the 
symbols mark the position of the snowball transition for several of the model planets with 24-hr days used in the sensitivity 
study discussed in Section 5.2. The orbital and planetary parameters for those planets are indicated above; for all other model 
planets P = 23.5°, 6 a = 0°, and f = 0.7. The black line shows where the (fixed-atmosphere) outer boundary of the habitable 
zone would be if it followed the (1 — e 2 ) -1 / 4 scaling that is sometimes assumed. This scaling has been used by previous studies 
to estimate the position of the habitable zone, and matches well with our modeled habitable zones for eccentricities < 0.65. 
Planets in more eccentric orbits appear to be habitable at greater semimajor axes than the simple scaling relationship would 
predict. For instance, the planet with 6 P = 90° and f = 0.7 is habitable out to 1.90 AU. Rotation rate has little effect on the 
position of the snowball transition for eccentricities between 0.35 and 0.65, but increasing rotation rate appears to extend the 
position of the habitable zone for highly eccentric planets. In low eccentricity orbits, the relationship between rotation rate and 
maximum habitable semimajor axis is not monotonic. The complex relationship may indicate that the semimajor axis of the 
snowball transition is maximized when the planet rotates slowly enough that sufficient heat is transported to the mid-latitudes 
to prevent the advance of global permanent ice coverage but quickly enough that a significant amount of heat is retained in the 
low- latitudes. 
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Fig. 9. — Temperature profile of a pseudo-Earth with rotation rate Q p = fi®, azimuthal obliquity a = 0, polar obliquity 
P = 90, and ocean fraction = 10% in an eccentric orbit with semimajor axis with a = 1.255 AU and eccentricity e = 0.20. 
Left: Habitability map of the planet with colors representing temperature as shown in the color bar. The planet is completely 
frozen for most of the year, but the southern hemisphere warms to just above 273 K during southern summer at periastron. 
Right: Temperature variations during the course of the year for the latitude bands at the equator, north pole, and south pole. 
Although the south pole is the only region of the planet that ever becomes habitable, it also experiences the most extreme 
temperature variations of any region on the planet. 



